Bounds for the refractive indices of metamaterials 
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1. Introduction 

During recent years several new types of artificial materials or metamaterials with 
sophisticated electromagnetic properties have been designed. The fabrication of custom 
structures with dimensions much smaller than the relevant wavelength has made it 
possible to tailor the effective electric permittivity e and the magnetic permeability //. 
For example, materials with thin metal wires simulate the response of a low-density 
plasma so that Ree may become negative in the microwave range pQ. Similarly, with 
the help of a split-ring structure, a strong magnetic resonance is achieved so that Re /i 
may be negative [2]. Passive media with Ree and Re/i simultaneously negative, first 
realized by Smith et al. are particularly interesting. Such materials are often referred 
to as left-handed, since, for negative e and the electric field, the magnetic field, and 
the wave vector form a left-handed set of vectors. As the Poynting vector and the wave 
vector point in opposite directions, the refraction at the boundary to a regular medium is 
negative. The concept of negative refraction, introduced by Veselago already in 1968 [I], 
has opened a new horizon of applications in electromagnetics and optics. In particular 
the possibility of manipulating the near-field may have considerable potential, enabling 
imaging with no limit on the resolution 

Materials with negative e and fi are necessarily dispersive |HJ E], and loss is 
unavoidable. Loss has serious consequences for the performance of certain components; 
for example, it has been shown that the resolution associated with the Veselago-Pendry 
lens is strongly dependent on the loss of the material IE] . Therefore, it is important 
to look for metamaterial designs with negative real part of the refractive index while 
the loss is low. In this paper, instead of performing a search in an infinite, complex 
design space, we will find ultimate, theoretical bounds based on causality. We will also 
find optimal e{uj) and h(uj) functions. For example, suppose our goal is refractive index 
close to -1 while loss is negligible in a limited bandwidth Q = [ui,u^i\. What is then 
the minimal variation of the refractive index in Q, given that the medium is passive? If 
we force the real part of the refractive index to be exactly -1 in Q, what will then be 
the minimal loss there? 

It is common to assume that the medium can be described by some specific model, 
such as, for example, single or multiple Lorentzian resonances. While this permits a 
straightforward analysis, it is not clear if a more general medium would give a more 
optimal response in some sense. We will therefore not use a spesific model, but rather 
assume only causality. 

2. Realizable electromagnetic parameters 

Any electromagnetic medium must be causal in the microscopic sense; the polarization 
and magnetization cannot precede the electric and magnetic fields, respectively. This 
means that e(oS) and fi(uj) obey the Kramers-Kronig relations. In terms of the 
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susceptibilities x = e ~l° r X = / i — 1> these relation can be written 

Imx = ftRex, (1) 
Re X = -Him X , (2) 

where TC denotes the Hilbert transform jHj. These conditions are equivalent to the fact 
that x is analytic in the upper half-plane (Imw > 0), and uniformly square integrable in 
the closed upper half-plane J. The susceptibilities are defined for negative frequencies 
by the symmetry relation 

X (-u) = (3) 

so that their inverse Fourier transforms are real. For passive media, in addition to 
©-© we have: 

Imx(w) >0 for uj > 0. (4) 

The losses, as given by the imaginary parts of the susceptibilities, can be vanishingly 
small; however they are always present unless we are considering vacuum [B]. 

Eqs. CE}-® imply that 1 + x is zero-free in the upper half-plane |B|. Thus the 
refractive index n = can always be chosen as an analytic function in the upper 

half-plane. With the additional choice that n — > +1 as u — > oo, n is determined 
uniquely, and it is easy to see that ©-(JIJ) hold for the substitution x — > n — 1. 

While any refractive index with positive imaginary part can be realized at a single 
frequency, the conditions CQ)-(@J) put serious limitations on what is possible to realize in 
a finite bandwidth. First we will investigate the possibility of designing materials with 
the real part of the refractive index less than unity. In particular we will analyze to 
what extent it is possible in a limited bandwidth to have a constant index below unity 
(or even below zero) while the loss is small. We set n — 1 = u + iv (or x — u + iv), where 
u and v are the real and imaginary parts of n — 1 (or x), respectively. To begin with, 
we set v(uj) = in the interval Q = [ux, cjz], (The case with a small imaginary part will 
be treated later.) By writing out the Hilbert transform and using (J3J), we find 

2 f"i v(u')u'du' 2 r 00 v(w')u}'du' 

+ ~ / — n ~7T ( 5 ) 



u{ uo 



7T JO LO 2 — Uj' 2 7T Jujo UJ 12 — LU 2 



for uj G f2. Note that both terms in (jHJ) are increasing functions of uj. Since the goal is 
a constant, negative u{ui) in Q, the second term should be as small as possible. In the 
limit where the second term is zero, we obtain 

. . 2 v (to')uj'du' u 2 - u 2 



and therefore 



,u 2 -u 2 



u{uj) — u{uj\) > \u(u>i)\ = — -, u G Q, (7) 



uJ 2 



| If the medium is conducting at zero frequency, the electric x is singular at lu = 0. Although \ is not 
square integrable in this case, similar relations as Q-© can be derived [Hj. 
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provided u(ui) is negative. In particular, the largest variation in the interval is 

u(uj 2 ) — u((Ji) > \u(ux)\(2A — A 2 ). (8) 

Here we have defined the normalized bandwidth A = (002 — ujx)juj%. These bounds are 
realistic in the sense that equality is obtained asymptotically when v(u) approaches a 
delta function in uo = + . In this limit u(uo) = u{uoi)uo 2 /uo 2 . 

It is interesting to estimate how much loss we must allow in the interval to wash 
out the variation (|7jl. Letting v(uo) approach a delta function in uo = + , and adding 



the positive function \u(ui)\J u 2 — oo 2 ^J oo 2 — uo\juo 2 for uo G VL correspond to a constant 
u(u>) = u(ui) in Q. Furthermore, it can be shown that this particular v corresponds 
to the minimal possible loss in the interval. The proof for this claim is given elsewhere 
|lUj . Thus the maximal value of the (minimal) loss in the interval satisfies 

^ > K^)!^^- = K^)I A + °( A2 )- ( 9 ) 

By a superposition of the optimal solutions associated with the bounds (J7j) and (JJJJ), 
we obtain a bound for the loss when a certain fraction 1 — a (0 < a < 1) of the variation 
(JJJ) remains: 

/ ,2 , ,2 

1 / \ 1 UO2 — 1 / \ 

^max > a w(o;i) — . (10) 

As an example, consider the case where the goal is refractive index close to -1 in an 
interval Q with A 1. In the limit of zero imaginary index, (JHJ) gives that the variation 
of the real index in the interval is larger than 4 A. It is interesting that the minimal 
variation is obtained approximately if the medium has sharp Lorentzian resonances for 
a low frequency. For example, let e{ui) = h{oj) = 1 + x{ u )i where 

TP, ) 2 

= 2 :° . r . (ii) 

(jJq — LO z — 10JL 



Here, F, ujq, and T are positive parameters. If the bandwidth T and center frequency ujq 
are much smaller than Ren(u;) ~ 1 — F(ujq/uj) 2 and Imn(o;) ~ Fuj 2 T/uj 3 for uo > u\. 
If we require Ren(wi) = —1, we obtain Imn(wi) ~ 2T/ui and Ren^) — Ren(ui) ~ 4A. 
When r/ui — > 0, this corresponds to the optimal refractive index function associated 
with the bound (JHJ. Furthermore, it is worth noting that if we want the real index 
variation to be zero in Q, the maximum imaginary part of the refractive index in Q must 
be larger than 2A. The required imaginary part in Q can roughly be approximated by 
weak resonances at {uj\ + uj 2 )/2, see Fig. 1. 

So far we have considered the case where the goal is a constant u(u) < in Q. If 
the goal is u(uj) > in Q, it is the last term in (jSJ) that comes to rescue. Inspired by 
the result (J2J), we may let u(u) approach a delta function at a frequency much larger 
than uj. Indeed, in the limit where this resonance frequency approaches infinity, the 
function u{uj) is constant and positive in VL while v{uS) is zero. Of course this limit is not 
realistic; in practice the resonance frequency is limited to, say u; max , where u; max > u>2- 
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Figure 1. The real (solid line) and imaginary (dashed line) refractive index associated 
with a Lorentzian resonance at w = ujq. The two figures represent the same functions 
but the scale is different. Also shown are the real part and imaginary part for the 
case where the real part is -1 in Q. In general this refractive index function can 
be found using the approach in |10j : however, when uj\ is much larger than the 
resonance frequency and bandwidth, the required v in ft is 2y 'lj 2 — uj 2 y/uj 2 — lu 2 /uj 2 . 
The parameters used are u>\ = 2wq, ^2 = 2.5woi T = O.Iuqi and F = 8. 



The associated bounds are easily deduced along the same lines as above. For example, 
(JHJ) becomes 

u) 2 — uj 2 

u(u 2 ) - u[uji) > tt(o;i)— y ^2- (12) 

Similarly, there may be a lower bound uj min , where < u; min < Ui, on the resonance 
frequency. The stricter inequalities in this case, corresponding to (fTjl- ffTUj) . can be found 
in a similar fashion. 

Eq. (jHJ) and (jl2j) have another interesting consequence. If the loss is zero in an 
infinitesimal bandwidth around u, we immediately find that the derivative du/du is 
bounded from below: 

du \2\u(u)\/u for u(u) < 0, r N 

- — > s 13 
did [0 for u{uj) > 0. 

For the case u{uj) > we have set c^Wx = oo. Note that also this bound is tight. A 
similar bound was obtained previously for e(u) and (jl(uj) [6J. Eqs. (|T3*jl should also be 
compared to the weaker bound dn/du > \u{uj)\/uj which was obtained recently [llj. 
While the latter bound means that the group velocities of transparent, passive media 
are bounded by c, (|T3|) implies the maximum group velocity c/ (2 — n) for n < 1 (and 
trivially c/n for n > 1). Here c is the vacuum light velocity. 
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When the loss in a bandwidth is at most v. 



max; 



(JT3J) becomes 



du 
du 




2\u(u)\ 4t>, 



TTUjA 



max 



) 



for u{uj) < 0, 



for u{uj) > 0, 



(14) 



to lowest order in A, for uj close to (loi +u> 2 )/2. In obtaining (fTl|) we have assumed that 
v(u) is approximately constant in Q, and calculated the corresponding contribution to 
du/du). A similar bound can be derived when v(uj) varies slowly in Q; for example, if 
v (lj) is the imaginary part of a Lorentzian with V = uj 2 — Ui, the inequality holds with 
the replacement 4/tt — > 2. Note that without an assumption on the variation of v(u), 
du/du) can take any value. 

A similar method as that leading to (jHJ can be used to find bounds for the variation 
of derivatives in Q, in the limit of no loss. For the first order derivative, the variation 
can be arbitrarily small to first order in A, for any positive du(ui)/du. For negative 
second order derivative D = d 2 u/du 2 (first order dispersion coefficient) we obtain 



3. Discussion and conclusion 

We have considered the set of realizable permittivities, permeabilities and refractive 
indices. For passive media we have used ©-© to prove ultimate bounds for the loss 
and variation of the real part of the permittivity, permeability, and refractive index. 

While the notation has indicated an isotropic medium, the bounds in this paper 
are valid for the effective index of the normal modes of anisotropic media as well. In 
this case the identification of associated e and fi tensors from the effective index is more 
complicated, but nevertheless feasible. More generally, the bounds are valid for the 
effective index of any electromagnetic mode that can be excited separately and causally, 
provided the effective index is independent of the longitudinal coordinate. 

On the basis of causality, it is clear that the susceptibilities of active media also 
satisfy (|T|) -([3"j). However, (J3J) is certainly not valid. Krem and Nudel'man have shown 
how to approximate a square integrable function in a finite bandwidth by a function 
satisfying CO)-© [El UHj- The approximation can be done with arbitrary precision; 
however, there is generally a trade-off between precision and the energy of x outside the 
interval [Hj . Once a valid susceptibility has been found, a possible refractive index can 
be found e.g. by setting n = e = fi=l + x- Hence, in principle, for active media n can 
approximate any square integrable function in a limited bandwidth. 
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